Current constraints on the dark energy equation of state 
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We combine complementary datasets from Cosmic Microwave Background (CMB) anisotropy 
measurements, high redshift supernovae (SN-Ia) observations and data from local cluster abun- 
dances and galaxy clustering (LSS) to constrain the dark energy equation of state parameter- 
ized by a constant pressure-to-density ratio wq. Under the assumption of flatness, we flnd 
wq < —0.85 at 68% c.l., providing no significant evidence for quintessential behaviour different 
from that of a cosmological constant. We then generalise our result to show that the constraints 
placed on a constant wq can be safely extended to dynamical theories. We consider a variety 
of quintessential dynamical models based on inverse power law, exponential and oscillatory 
scaling potentials. We find that SNla observations are 'numbed' to dynamical shifts in the 
equation of state, making the prospect of reconstructing w{z), a challenging one indeed. 



Introduction. The discovery that the universe's evo- 
lution may be dominated by an effective cosmological 
constant , is one of the most remarkable cosmolog- 
ical findings of recent years. An exceptional opportu- 
nity is now opening up to decipher the nature of dark 
matter |^ , to test the veracity of theories and recon- 
struct the dark matter's profile using a wide variety 
of observations over a broad redshift range. 

One matter candidate that could possibly explain 
the observations is a dynamical scalar "quintessence" 
field. One of the strong aspects of quintessence theo- 
ries is that they go some way to explaining the fine- 
tuning problem, why the energy density producing the 
acceleration is ~ 10~^^''Mp;. A vast range of "tracker" 
(see for example |^,^) and "scaling" (for example ||]- 
11]) quintessence models exist which approach attrac- 
tor solutions, giving the required energy density, inde- 
pendent of initial conditions. The common character- 
istic of quintessence models is that their equations of 
state, iuq = p/ P, vary with time whilst a cosmological 
constant remains fixed at iuq^a = — 1. Observation- 
ally distinguishing a time variation in the equation of 
state or finding wq different from —1 will therefore be 
a success for the quintessential scenario. 

In this paper we will combine the latest observa- 
tions of the Cosmic Microwave Background (CMB) 
anisotropies provided by the Boomerang ||], DASI 
ficjl and Maxima experiments and the informa- 
tion from Large Scale Structure (LSS) with the lumi- 
nosity distance of high redshift supernovae (SN-Ia) to 
put constraints on the dark energy equation of state 
parameterised by a redshift independent quintessence- 
field pressure-to-density ratio wg.We will also make 
use of the Hubble Space Telescope (HST) constraint 
on the Hubble parameter h = 0.72 ± 0.08 We 
will then also consider whether one can feasibly ex- 
tract information about the time variation of w from 
observations. 

The importance of combining different data sets in 
order to obtain reliable constraints on wq has been 
stressed by many authors (see e.g. jl3|, [|l5|, fl^), 
since each dataset suffers from degeneracies between 
the various cosmological parameters and wq . Even 



if one restricts consideration to flat universes and to 
a value of wq constant with time the SN-Ia luminos- 
ity distance and position of the first CMB peak are 
highly degenerate in wq and Qq, the energy density 
in quintessence. 

The paper is therefore structured as follows: in sec- 
tion II and III we will present the CMB, SN-Ia and 
LSS data used in the analysis. In section IV we will 
present the results of our analysis. We will consider 
the implications for a dynamical wq in section V. Sec- 
tion VI will be devoted to the discussion of the result 
and the conclusions. 

Constraints from CMB. The effects of quintessence 
on the angular power spectrum of the CMB 
anisotropies are several [p5 17 , l3]. In the class of mod- 



els we are considering, however, with a negligible value 
of Qq in the early universe 

in order to satisfy the BBN bound ( ||l^), the effects 
can be reduced to just two . 

Firstly, since the inclusion of quintessence changes 
the overall content of matter and energy, the angu- 
lar diameter distance of the acoustic horizon size at 
recombination will be altered. 

In flat models (i.e. where the energy density in mat- 
ter is equal to Om = 1 — ^q), this creates a shift in 
the peaks positions of the angular spectrum as 
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It is important to note that the effect is completely 
degenerate in the interplay between wq and Qq. Fur- 
thermore, it does not qualitatively add any new fea- 
tures additional to those produced by the presence of 
a cosmological constant and it is not highly sen- 
sitive to further time dependencies of wq . 

Secondly, the time-varying Newtonian potential af- 
ter decoupling will produce anisotropies at large angu- 
lar scales through the Integrated Sachs- Wolfe (ISW) 
effect. The curve in the CMB angular spectrum on 
large angular scales depends not only on the value of 
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FIG. 1. CMB power spectra and the angular diameter 
distance degeneracy. The models are computed assuming 
flatness, flfc = 1 — SIm — = 0). The Integrated Sachs 
Wolfe effect on large angular scale slightly breaks the de- 
generacy. The degeneracy can be broken with a strong 
prior on h, in this paper we use the results from the HST. 

wq but also its variation with redshift. However, this 
effect will be difficult to disentangle from the same ef- 
fect generated by a cosmological constant, especially 
in view of the affect of cosmic variance and/or gravity 
waves on the large scale anisotropics. 

In order to emphasize the importance of degenera- 
cies between all these parameters while analyzing the 
CMB data, we plot in Figure 1 some degenerate spec- 
tra, obtained keeping the physical density in matter 
r^M^^, the physical density in baryons ftbh^ and TZ 
fixed. As we can see from the plot, models degenerate 
in Wq can be constructed. However, as we will utilise 
in the next sections, the combination of the different 
datasets can break the mentioned degeneracies. 

To constrain wq from CMB, we perform a like- 
lihood analysis comparing the recent CMB observa- 
tions with a set of models with cosmological param- 
eters sampled as follows: 0.1 < < 1-0, —1.0 < 
WQ < -0.55, 0.015 < Qb < 0.20; < Qq < 0.9 
and 0.45 < h < 0.95. We vary the spectral index of 
the primordial density perturbations within the range 
Us = 0.60, 1.40, we allow for a possible reionization 
of the intergalactic medium by varying the CMB pho- 
tons optical depth in the range 0.0 < tq < 0.4 and 
we re-scale the fluctuation amplitude by a pre-factor 
Cio, in units of Ciif^^. We also restrict our analysis 
to flat models, fltot = 1, and we add a conservative 
external prior on the age of the universe to > 10 Gyrs 
(see e.g. [p4[). 

In order to speed-up the computation time of the 
theoretical models for different wq we make use of a 
/c-splitting technique |20|. Basically a wq = — 1 and 
WQ > — 1 model are calculated in two different ways. 
For low £ the 2 models are computed in the ordinary 



FIG. 2. Contours of constant R (CMB) and SN - la 
luminosity distance in the wq-Q.m plane. The degeneracy 
between the two distance measures can be broken by com- 
bining the two sets of complementary information. The 
luminosity distance is chosen to be equal to d; at z = 1 for 
a fiducial model with = 0.7, Q,m = 0.3,/i = 0.65. (We 
note that as JIq = 1 — Q.m goes to zero the dependence of 
R and upon wq also become zero, as there is no dark 
energy present.) 

way by solving the Boltzmann equation, in order to 
properly take in to account the ISW effect. For the 
high I just a ffat, wq — —1 model is calculated. This 
Wq = — 1 model is then shifted using the expression 
for the angular diameter distance in equation (1) to 
obtain the wq > — 1 models. 

The theoretical models are computed using a mod- 
ified version of the publicly available CMBFAST pro- 
gram ( p^) and are compared with the recent 
BOOMERanG-98, DASI and MAXIMA-1 results. 
The power spectra from these experiments were esti- 
mated in 19, 9 and 13 bins respectively, spanning the 
range 25 < ^ < 1100. We approximate the experimen- 
tal signal inside the bin to be a Gaussian variable, 
and we compute the corresponding theoretical value 
by convolving the spectra computed by CMB- 
FAST with the respective window functions. When 
the window functions are not available, as in the case 
of Boomerang-98, we use top-hat window functions. 
The likelihood for a given cosmological model is then 
defined by -2\nC = [C^^ - C%'')MBB'iC'j^> - Cf,) 
where C|j'' (C^^) is the theoretical (experimental) 
band power and Mbb' is the Gaussian curvature of 
the likelihood matrix at the peak. We consider 10%, 
4% and 4% Gaussian distributed calibration errors (in 
^iK)ior the BOOMERanG-98, DASI and MAXIMA- 
1 experiments respectively and we take in to account 
for the beam error in BOOMERanG-98 by analytic 
marginalization as in |Q . We also include the COBE 
data using Lloyd Knox's RADPack packages. 

Constraints from Supernovae. Evidence that the 
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TABLE I. Constraints on wq and Q,m = 1 — Qq using 
different priors and datasets. We always assume flatness 
and to > 10 Gyr. Tfie Icr limits are found from the 16% 
and 84% integrals of the marginalized likelihood. The HST 
prior is h = 0.72 ±0.08 while for the BBN prior we use the 
conservative bound »b/i^ = 0.020 ± 0.005. 



CMB+HST WQ < -0.62 

0.15 < riM < 0.45 
CMB+HST+BBN -0.95 < wq < -0.62 

0.15 < < 0.42 
CMB+HST+SN-Ia -0.94 < wq < -0.74 

0.16 < < 0.34 
CMB+HST+SN-Ia+LSS wq < -0.85 

0.28 < riM < 0.43 



universe's expansion rate was accelerating was first 
provided by two groups, the SCP and High-Z Search 
Teani( |l]) using type la supernovae (SN-Ia) to probe 
the nearby expansion dynamics. SN-Ia are good stan- 
dard candles, as they exhibit a strong phenomenologi- 
cal correlation between the decline rate and peak mag- 
nitude of the luminosity. The observed apparent bolo- 
metric luminosity is related to the luminosity distance, 
measured in Mpc, hy ms = M +5logdL{z)+25. where 
M is the absolute bolometric magnitude. The lumi- 
nosity distance is sensitive to the cosmological evolu- 
tion through an integral dependence on the Hubble 
factor di = {1 + z) {dz' / H{z' ,V,Q,Wq) and there- 
fore can be used to constrain the scalar equation 
of state. We evaluate the dark energy fl/w likeli- 
hoods assuming a constant equation of state, such 
that H{z) = poE.if^«(l + The predicted 

rrief f is then calculated by calibration with low-z su- 
pernovae observations p5| ] where the Hubble relation 
di ~ Hqcz is obeyed. We calculate the likehhood, £, 
using the relation C = £o exp(— X^(^j ''^q)/2) where 
Co is an arbitrary normalisation and is evaluated 
using the observations of the SCP group, marginalis- 
ing over Hq . As can be seen in Figure || there is an 
inherent degeneracy in the luminosity distance in the 
^Im/wq plane; one can see that little can be found out 
about the equation of state from luminosity distance 
data alone. However, the degeneracies of CMB and 
SNla data complement one another so that together 
they offer a more powerful approach for constraining 

WQ- 

Results. Table I shows the 1-a constraints on wq 
for different combinations of priors, obtained after 
marginalizing over all remaining nuisance parameters. 
The analysis is restricted to flat universes. One can see 
that Wq is poorly constrained from CMB data alone. 
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FIG. 3. The likelihood contours in the {Hm, wq) plane, 
with the remaining parameters taking their best-fitting 
values for the joint CMB-f SN-Ia-|-LSS analysis described 
in the text. The contours correspond to 0.32, 0.05 and 0.01 
of the peak value of the likelihood, which are the 68%, 95% 
and 99% confidence levels respectively. 

even when the HST strong prior on the Hubble param- 
eter, h = 0.72 ± 0.08, is assumed. Adding a Big Bang 
Nucleosynthesis prior, ilbh^ — 0.020±0.005, has small 
effect on the CMB-hHST result. Adding SN-Ia breaks 
the CMB Qq — WQ degeneracy and improves the up- 
per limit on wq, with wq < —0.74. Finally, includ- 
ing information from local cluster abundances through 
(Tg = {0.55±0.1)Q~^j '^, where as is the rms mass fluc- 
tuation in spheres of 8h~^ Mpc, further breaks the 
quintessential-degeneracy, giving wq < —0.85 at 1- 
a. Also reported in Table I, are the constraints on 
flAi- As we can see, the combined data suggests the 
presence of dark energy with high significance, even in 
the case CMB-I-HST. It is interesting to project our 
likelihood in the fig — wq plane. Proceeding as in 
p^ , we attribute a likelihood to a point in the (JIm, 
Wq) plane by finding the remaining parameters that 
maximise it. We then define our 68%, 95% and 99% 
contours to be where the likehhood falls to 0.32, 0.05 
and 0.01 of its peak value, as would be the case for 
a two dimensional multivariate Gaussian. In Figure 3 
we plot likelihood contours in the {^Im, wq) plane for 
the joint analyses of CMB-^SN-Ia-hHST-t-LSS data to- 
gether with the contours from the SN-Ia dataset only. 
As we can see, the combination of the data breaks the 
luminosity distance degeneracy. 

Constraining dynamical models. In the previous 
section we obtained bounds on the equation of state 
parameter wq by assuming it is a constant, indepen- 
dent of the redshift. However, in quintessential models 
the equation of state can vary with time. It is there- 
fore useful to discuss how well our constraints on wq 
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apply to less trivial models. There are a wide vari- 
ety of quintessential models; we illustrate our analy- 
sis using representatives of three of the most general 
classes of model: the inverse power law,y((/)) — Vb/f/'^ 
1^ with p=l, an exponential scaling potential with 
a feature, V{<^) = V^e'^'^^A -h ((/i - 0o)^) H], with 
A = 10, ^ = 0.008, 00 — 25.8, and an oscillatory scal- 
ing potential,l/(0) = V^e-^'^il + Asin(y(^)) ^ with 
A = 4, A = 0.98,1^ = 0.51. In each case Vb is chosen 
so that JIq = 0.7 and is 65. The particular time 
dependent characteristics of each are shown in Figure 



show that even a substantial deviation from We// at 
late times produces a small change in pq. The ef- 
fect on ms is doubly 'numbed', firstly because of the 
smoothing by the integral relation with H in d;, see 
Q for a previous discussion of this, and secondly be- 
cause rriB ^ logdi. As a result the bolometric magni- 
tude taken from the SNla data is highly insensitive to 
variations in the equation of state. This doesn't bode 
well if we are to try and reconstruct the time vary- 
ing equation of state from observations. It looks more 
likely that Weff will be a more tangible observable. 




FIG. 4. The variation of Qq and wq with redshift for 
the three models described in section V. The power law po- 
tential (full line) WQ shows a steady small variation, the 
exponential feature potential, (short dashed line) acts re- 
markably like a cosmological constant at late times, after a 
deviation from behaving like normal matter with wq ~ 0, 
whilst the oscillating potential (dot-dash) shows a contin- 
ual variation in wq. 

If we are to constrain dynamical models, we need 
to understand how well the effect of a time varying 
can be modeled by a constant Weff- In |28| it 
was shown that in models in which the dark energy 
component is negligible at last scattering, the CMB 
and matter power spectra are well approximated by 
Weff — J dartQ{a)wQ{a). This is demonstrated in fig- 
ure |5| in the case of the exponential potential with a 
feature, using Weff = —0.993. In we showed that 
if the dark energy component is a significant propor- 
tion, as can be seen in scaling quintessence models, 
the dark energy component can be modelled as an 
additional contribution to the effective number of rel- 
ativistic degrees of freedom. We restrict ourselves to 
former case in which Qfj,{AIeV) is negligible. In ||l3| 
it is noted that although Weff is a good measure for 
modelling CMB spectra it may not be such a good 
measure when considering di . We investigate whether 
this is actually the case using the three models out- 
lined above, a power law, exponential scaling and scal- 
ing oscillatory potential, as test cases. In Figure o we 
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FIG. 5. Comparison of CMB temperature power spec- 
tra for the dynamical quintessence model with an expo- 
nential potential with a feature described in section V, 
and a model with constant wq = We// for the dynamical 
model. In both cases with Qq — 0.7 and Hq = 65. One 
can see that the constant model is a remarkably good ap- 
proximation to the dynamical model, despite the equation 
of state of the dynamical model varying significantly from 
the effective value from recombination until nowadays. 



Conclusions. In this paper we have provided new 
constraints on the dark energy equation of state pa- 
rameter Wq by combining different cosmological data. 
The new CMB results provided by Boomerang and 
DASI improve the constraints from previous and sim- 
ilar analysis (see e.g., |29|), with wq < —0.85 
at 68% c.l. {wq < -0.76 at 95% c.L). We have also 



demonstrated how the combination of CMB data with 
other datasets is crucial in order to break the ^Iq — wq 
degeneracy. The constraints from each single datasets 
are, as expected, quite broad but compatible between 
each other, providing an important consistency test. 
When comparison is possible (i.e. restricting to sim- 
ilar priors and datasets), our analysis is compatible 
with other recent analysis on wq ( pH). Our final 
result is perfectly in agreement with the wq = — 1 
cosmological constant case and gives no support to a 
quintessential field scenario with wq > — 1. A frus- 
trated network of domain walls or a purely exponen- 
tial scaling field are excluded at high significance. In 
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FIG. 6. The imphcation of a deviation 
AwQ = 1 — w^ff/wQ on ADl = 1 ~ DL{weff)/DL{wQ) 
and Arrifl — mB(wQ) — mB{weff) for the three mod- 
els in section V. The traclcer potential (full line) has 
Weff ~ —0.727, the exponential potential with feature 
(dashed line) has Weff ~ —0.993 and the oscillatory po- 
tential (dot-dash) has Weff = —0.529. Notice that Dl and 
TUB are in turn both desensitized to any variation in wq. 



addition a number of quintessential models are highly 
disfavored, power law potentials with p > 1 and the 
oscillatory potential discussed in this paper, to name 
a few. 

It will be the duty of higher redshift datasets, for 
example from clustering observations |Q to point to 
a variation in w that might place quintessence in a 
more favorable light. 

The result obtained here, however, could be plagued 
by some of the theoretical assumptions we made. The 
CMB and LSS constraints can be weakened by the in- 
clusion of a background of gravity waves or of isocur- 
vature perturbations or by adding features in the pri- 
mordial perturbation spectra. These modifications 
are not expected in the most basic and simplified in- 
flationary scenario but they are still compatible with 
the present data. The SN-Ia result has been obtained 
under the assumption of a constant- with-time wq . We 
have shown that in general Weff is a rather good 
approximation for dynamical quintessential mode Is 
since the luminosity distance depends on wq through 
a multiple integral that smears its redshift depen- 
dence. As we show in the previous section, our re- 
sult is therefore valid for a wide class of quintessential 
models. This 'numbing' of sensitivity to wq implies 
that maybe an effective equation of state is the most 
tangible parameter able to be extracted from super- 
novae. However with the promise of large data sets 
from Planck and SNAP satellites, opportunities may 
yet still be open to reconstruct a time varying equa- 
tion of state m&i . 
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